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Abstract 

We propose a general theoretical scheme for the investigation of light emitted from nano- and 
micrometric structures of arbitrary shape and composition. More specifically, the proposed fully 
three-dimensional approach allows to derive the light-intensity distributions around the emitting 
structures and their modifications in the presence of nearby scattering objects. Our analysis allows 
to better identify the non-trivial relationship between near-field images and fluorescent objects. 

PACS numbers: 78.67.-n, 42.50.Ct, 07.79.Lh 
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Recent continuous progress in scanning near-field microscopy together with the develop- 
ment of adequate nanofabrication techniques has enhanced our insight into the field distri- 
butions in the proximity of nano- and microstructured materials. This insight is of great 
relevance for the design of future optical circuitry able to process light with the versatility of 
electronic chips. These developments have stimulated more refined theoretical approaches 
as well as novel simulation strategies for the analysis of light propagation in complex struc- 
tures and have renewed the interest in the classical theory of light scattering [I], ^ . Indeed, 
theoretical predictions are essential for the design of photonic structures able to control the 
propagation of light. As a matter of fact, most of the theoretical investigations performed 
so far focus on the scattering of incident light illuminating structures of non-trivial shapes 
|l|, [|, or study light modes in photonic crystals and microcavities ||. However, optical 
circuitry, besides passive elements controlling/ manipulating the flow of photons, require ac- 
tive components able to emit and/or amplify light. This, in turn, opens relevant questions 
concerning the emission patterns of active mesoscopic systems: Which is the actual field 
distribution around emitting structures? How does the interaction with nearby scattering 
objects modify the field distribution? What kind of relationship between near-field images 
and fluorescent objects holds? 

Aim of this Letter is to provide a general theoretical framework for the evaluation of 
the field distributions in the proximity of three-dimensional (3D) mesoscopic fluorescent 
objects in the presence of nearby scattering structures. The proposed quantum theory of 
light emission can be applied to a great variety of nanostructured optically active materials 
as mesoscopic dielectric objects uniformly doped with optically active molecules (e.g., dye 
molecules) or embedding semiconductor dots or layers able to emit light when appropriately 
excited. 

As a starting point, let us consider the key quantity we want to investigate, i.e., the 
spectrally-resolved energy density X(r, u>) corresponding to the electric field at point r; this 
can be defined as: 



where E + (r, to) is the electric- field operator corresponding to the positive frequency u (it can 
be expanded in terms of photon destruction operators) and E~(r,c<j) is its Hermitian con- 
jugate. Our theoretical approach is based on the Green's tensor technique in the frequency 
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domain |], 0- The electric field operator at a given positive frequency can be regarded as 
the sum of scattering and emission contributions ||: 

E+(r,^)=E+(r, W ) + E+(r,^) . (2) 

The first term describes light coming from free-space and scattered by the material system, 
i.e., 

E+(r, u) = E+(r, u) - k 2 J G(r, r', u) X (r', w)E+(r', u)dr' , (3) 

where k = u/c, the free-space electric-field operator Eq describes input light, and x = e — \ 
is the susceptibility function of the material system (that we have assumed to be a local and 
scalar function to avoid complications), and G is the Green tensor obeying the following 
equation 

-V x V x G(r, r , u) + k 2 £{r, w)G(r, r, u) = I5(r, r') , (4) 

where I is the unit dyadic. We stress that nonlocality and/or anisotropy can be easily 
implemented ||. The second term in Eq. (|2]) describes light emitted by the material system 
itself, i.e., 

E+(r,w) = -iw/io J G(r,r',^)j(r',^)rfr', (5) 

where the integration is performed over the volume of the scattering system and /xq denotes 
the magnetic permeability of vacuum. Here, the quantum noise operators j are the sources of 
spontaneous light emission. These zero-mean operators can be derived from the Heisenberg- 
Langevin equations for the material system || and are present only if the susceptibility 
tensor describing the material system has an imaginary part; they are a direct consequence 
of the fluctuation-dissipation theorem and obey the following commutation rules: || 

ji(r,u),jif(r',uj) =0, (6) 
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Ix^^l^^r-rO^-o;'), (7) 



where x 1 indicates the imaginary part of x- By inserting Eq. (0) into Eq. ([I]), it is easy to 

verify that the intensity of the emitted light is directly related to the correlation function of 

quantum noise currents: (j^ • j^. In order to calculate it explicitly we adopt a simple 2- level 

molecular model to describe the dielectric function of the scattering system, 

./v x d 2 N a(r) - N b (r) . . 

e h oo — uj + 27 



FIG. 1: (a) Field-intensity distribution around a pair of nanoscopic emitting quavers with uniform 
composition and population densities; (b) as (a) but with the elliptic pad on the right of the 
structure replaced with a silver pad. The two calculations have been made at A = 600 nm and on 
the observation plane located 30 nm above the structure. 

where Xs is the real susceptibility of the scattering system in the absence of active molecules, 
d is the transition dipole moment, M gives the density of active molecules, and N a ^) indicates 
the population densities of the upper (a) and lower (b) levels. This simple model can be easily 
extended to include non-random molecular orientations, more energy levels and molecular 
species. Also more complex models describing, e.g., semiconductor dielectric functions, 
can be adopted. By using this simple model we obtain, 

(j/(r,^)jHrV)) = —^7( I {T,Lj)N{r)S u ,6{T-i f )S{LJ-Ly), (9) 

with the population factor N(r) = N a (r)/(N b (r) — N a (r)). In the following we consider the 
case of absence of input photons (CEq(t,u) ■ Eq (r,u/)) = 0). By combining Eqs. (Q) and 
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FIG. 2: (a) Field-intensity distribution around a pair of nanoscopic emitting quavers with uniform 
composition and population densities; (b) as (a) but with the elliptic pad on the right of the 
structure replaced with a silver pad. The two calculations have been made at A = 600 nm and on 
the observation plane located 30 nm above the structure. 



we obtain I(r,u) = ^Z;(r,co>) with 

hk 4 r 

X t (r,u;) = —J2 / iV(rV(r>) |G^(r,rV)| 2 dr'. 
ii J 



(10) 



We observe that population densities in principle are affected by light emission, thus Eq. 
(|ITi|) should be solved self-consistently together with the rate equations for the populations. 
However, in many cases (far from the laser threshold), population densities are poorly af- 
fected by light emission and propagation, being fixed by external pumping. We also observe 
that Eq. ( |10|) has a structure analogous to the corresponding equation describing emission 
from thermal sources [|7| || . 

To assess the power and versability of the proposed scheme, we investigate the field- 
intensity distribution around a fluorescent low-simmetry nano-structured system. We con- 
sider first a pair of 3D emitting quavers with uniform composition (relative dielectric constant 
e = 8.9 + 0.5i) and population densities. The size of the whole structure is 450 x 380 nm 
and its height is 30 nm. Fig. 1 (a) displays I(x, y, z, u) in the proximity of the emitting 





FIG. 3: The field intensity I(r,u) (a) and partial intensities 2/(r,u;) for the structure with the 
metallic pad; (b) l x , (c) l y , (d) I z . 

object in an observation plane located at z = 30 nm. The Green tensor is computed by 
solving numerically a discretization of the Dyson equation JlJ . Calculations have been made 
at A = 2nc/u = 600 nm. The emission pattern reproduces almost perfectly the object 
shape (see also Fig. 2a where the structure is outlined). The light intensity is highest at 
the structure edges and nearby the smaller structures. Fig. lb shows the effect of replacing 
the elliptic pad on the right of the structure with a silver pad. The emission pattern results 
to be significantly altered. In particular a strong enhancement of the field intensity nearby 
the edge of the silver pad is clearly observed while intensity is strongly depleted towards 
the pad center. This behaviour is typical of surface plasmons ||, in this case excited by the 
incoherent light originating from the nearby fluorescent structure. 

The vector character of the emitted light is shown in Fig. 2. It displays the partial 
intensities 2j(r, uS) (i.e. the intensities that can be detected by a probe able to select light 
polarization) for the structure with the metallic pad. In spite of the isotropic character of 
the source currents, the three intensity distributions strongly differ one from each other. 
We observe that only X z reproduces quite well the object shape (except the metallic pad 
of course). Fig. 3 displays l(x,y,z,u) for (a) z = 30, (b) z = 60, (c) z = 90, and (d) 



6 



z = 150 nm above the object. It shows how the relationship between near-field images and 
fluorescent objects is progressively lost at increasing observation distances. 

We have proposed a general theoretical scheme for the investigation of light emitted from 
nano- and micrometric structures of arbitrary shape and composition. The presented numer- 
ical results provide useful guidelines for the interpretation of near-field photoluminescence 
spectroscopy/microscopy measurements. Moreover this theoretical and numerical scheme 
can be applied to the design of photonic systems with active elements. 

We thank Omar Di Stefano for stimulating discussions. 
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FIGURE CAPTIONS 



Fig. 1 

(a) Field-intensity distribution around a pair of nanoscopic emitting quavers with uniform 
composition and population densities; (b) as (a) but with the elliptic pad on the right of 
the structure replaced with a silver pad. The two calculations have been made at A = 600 
nm and on the observation plane located 30 nm above the structure. 



Fig. 2 

The field intensity I(r,u) (a) and partial intensities Xi(r,u) for the structure with the 
metallic pad; (b) I x , (c) X y , (d) l z . 



Fig. 3 

The field intensity X calculated on observation planes at increasing distance from the 
structure; (a) z— 30 nm, (b) z— 60 nm, (c) z— 90 nm, (d) z— 150 nm. 
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